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Abstract. It is a consequence of QCD's spontaneously broken chiral symmetry that 
the response of the neutron to long wavelength probes will be driven by excitation 
of its pion cloud. In this article I discuss manifestations of this in electromagnetic 
form factors and Compton scattering. Particular attention is given to the interplay 
between this physics and shorter-distance components of the neutron's electromagnetic 
structure. 
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1. Introduction 

In the limit of massless up and down quarks the QCD Lagrangian has an SU{2)l x 
SU{2)ji symmetry. We know from lattice simulations that the vacuum of QCD in this 
limit does not share this symmetry of the Lagrangian, with the axial combination of 
S'[/(2)'s being spontaneously broken [Ij. This confirms the key elements of the picture 
postulated by Nambu in talks and papers of 1960-1, work for which he was awarded 
the 2008 Nobel Prize in Physics [2l|3]. An immediate consequence of the spontaneous 
breaking of chiral symmetry is that there will be three Goldstone bosons: the tt^, tt^, 
and 7T~ . Since the up and down quarks are not exactly massless these pions are also not 
massless, but, to the extent that rriu and are small compared to Aqcd^ the pions 
will be much lighter than all other hadronic degrees of freedom. 

This picture is systematized in an effective field theory (EFT) called 'chiral 
perturbation theory' (xPT) [4[ [5l |6l E]. Chiral perturbation theory encodes the 
SU{2)l X SU{2)pi symmetry of QCD, and the pattern of its breaking, in a Lagrangian 
containing nucleon and pion fields. It includes all possible interactions consistent with 
the symmetry. These are organized in a derivative expansion, together with an expansion 
in powers of rUq for the symmetry-breaking operators. The EFT can be used to calculate 
both tree and loop diagrams, which leads to a perturbation expansion for scattering 
amplitudes in powers of 

p^P^_ (1) 

A 

Here ^^fir ^ 1200 MeV (with f^^ ^ 93 MeV the pion-decay constant) is the scale 

that emerges from loop calculations [8j as controlling the convergence of the expansion 
for meson-meson interactions. If we consider threshold processes, and so p = 0, this 
is an expansion in powers of m^^ and so can be considered an expansion around the 
chiral limit of QCD, where tUt^ = 0. In the meson sector the theory is mature, and 
calculations are routinely performed to two loops. For nn scattering the results are 
extremely accurate, reaching the sub-1% level. (For a recent review see [9].) 

Adding baryons to the theory as static sources facilitates a power counting in 
which each diagram has a definite order in the P expansion [lOj. The effects of baryon 
recoil constitute corrections to the static-source picture, and they can be systematically 
included in the Lagrangian via an expansion in powers of p/M pHl [TT j. The result is 
known as "heavy-baryon chiral perturbation theory" (HB^PT). HB^PT for nucleons 
and pions is summarized in Section [2} There has been a recent fiurry of activity 
devoted to baryon xPT calculations that do not invoke this p/M expansion [121 fT3] . 
but here I will be mainly using the heavy-baryon variant of the theory to describe 
neutron structure. Electromagnetic interactions are easily included in this theory, via 
minimal substitution and the addition of terms that are gauge invariant by themselves 
(e.g. an operator corresponding to magnetic fields coupling to the neutron anomalous 
magnetic moment). In the usual xPT counting the electron charge e then counts as one 
power of P. For comprehensive reviews of xPT in the baryon sector (both heavy-baryon 
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and "relativistic" formulations) see Refs. p!4l [T5l [T6l [T7] . 

The standard xPT Lagrangian includes only nucleons and pions, and so its 
breakdown scale is not 4:7Tf^. Instead it is set by the lightest hadronic degree of freedom 
not explicitly included in the theory, which in most reactions is the Delta(1232). So, 
unless the xPT Lagrangian is extended to include Delta degrees of freedom, the scale 
in the denominator in Eq. ([T]) is not A^, but is instead Ma — ^ 290 MeV. While 
there has been much progress in this direction [TOl [HI [IHl EOj EH [22] , for most of this 
article I will regard the Delta-nucleon mass difference as a high-energy scale. 

My focus here is thus on the chiral structure of the neutron, as revealed by 
electromagnetic probes. I therefore report on computations of electron and photon 
scattering, with the main emphasis on results for these reactions obtained in HB^PT. 
I will be particularly concerned with the interplay between pion-cloud mechanisms and 
the role of higher-energy excitations. The ability to vary photon energy, cj, in the case of 
Compton scattering, and three-momentum transfer, q, in the case of electron scattering, 
allows us to map out this interplay over a significant kinematic domain. It also provides 
interesting data on the breakdown of xPT, since comparison of its predictions with 
experimental data over a range of u and |q| can reveal where the theory is failing, 
and what additional degrees of freedom are needed in order to extend its radius of 
convergence. 

But, practically speaking, these aspects of neutron structure cannot be examined 
separately from the chiral structure of nuclei. Neutrons cannot be collected in sufficient 
numbers to make electromagnetic scattering reactions on free neutrons an experimental 
reality. Hence everything we know about neutron electromagnetic structure is inferred 
from more complex processes, e.g. experiments on light nuclei. In order to extract 
information on neutron structure from such experiments an understanding of the multi- 
nucleon effects that are present in the data is also necessary. 

Developments over the past 15 years have also allowed xPT to provide a systematic 
approach to multi-nucleon problems. Weinberg's seminal papers on xPT expansions for 
the nuclear potential [23l [2l] and the operators governing the interaction of external 
probes with the nucleus [25] facilitate a parallel and consistent xPT treatment of neutron 
structure and electromagnetic reactions in two- (and three- and . . . ) nucleon systems. 
In Section [3] I will give a brief explanation of how xPT enables the consistent treatment 
of the operators governing processes in single- and multi-nucleon systems. Electron- 
neutron scattering and jn jn will then be discussed in Sections [4] and [5] respectively. 
In this article I only have space to discuss those two processes in detail, so in the 
final section. Section |6| I briefiy mention some other reactions that probe the same 
physics: virtual Compton scattering, as well as pion photo- and electro-production near 
its threshold. I also offer some concluding thoughts on the neutron's chiral structure. 
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2. Chiral perturbation theory for nucleons and pions 

In this section I provide a brief summary of the aspects of chiral perturbation theory 
(xPT) that are relevant for the subsequent presentation. The discussion here follows 
the notation and conventions of Ref. [14]. In particular, I will use the 'heavy-baryon' 
formulation of xPT, in which a p/M expansion is made alongside the expansion in 
powers ofp/A^. This expansion is not essential to the efficacy of xPT, but since M ^ 
it does not have a deleterious effect on its convergence — apart from a few specific cases 
that involve proximity to kinematic thresholds. 

XPT is an effective field theory of QCD, that encodes the symmetries of the 
underlying theory which are pertinent to the regime where energies are of the order 
of Goldstone boson (here, pion) masses. As such the chiral Lagrangian is formulated in 
terms of objects with straightforward transformation properties under the chiral group 
SU (2) L X SU For details on the transformation properties of the building blocks of 
the chiral Lagrangian the reader is referred to Refs. [HI [26]. Here we choose a particular 
representation for the pion field, and summarize the building blocks as: 

u = exp (^^y^ j ; U = u^] (2) 

where tt is the pion iso-triplet. From these we construct the chiral covariant derivative 
of the nucleon field, L^^, and the axial- vector object with one derivative, Uj^: 

D^^ = d^^ + r^T/;; (3) 

^ L^dfj^u + udfj^v) — ieAfj^ [y) Zu + uZv)^ ; (4) 

u^d^u — ud^v) — leA^ (y^Zu — uZv)^ , (5) 

with 2^ = |(1 + ra). Here we have specialized to the external-field case of interest for 
this paper: A^ is an (external) photon field and there is no external axial field. Note 
that the vector field F^ contains only even powers of the pion field, and contains only 
odd powers of it. The nucleon field is the usual iso-double. Both components of ij) 
are Dirac fields. 

The leading-order Lagrangian contains all possible terms with zero and one 
derivative and is: 

= ^ (. ^ - M + ^7^«p75) ^. (6) 

Strictly speaking M and qa here are the mass and axial coupling of the nucleon in the 
chiral limit, although the distinction between their values at zero quark mass and the 
physical ones will not be important for any of the discussions here. This Lagrangian, 
supplemented by the standard leading-order meson-meson one: 

= ^-^[Ti{D^U^D^U) - 2BTt{M{U + U^))] (7) 

with Df^U = d^U — ieAfj^l^Ts^U]^ is all one needs to calculate the leading-order loop 
effects in chiral perturbation theory. 
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Such calculations are rendered more straightforward, and more transparent, if one 
removes the large energy scale associated with the nucleon mass. This energy is not 
available for particle production, because of baryon-number conservation. Hence it 
plays no role in the single- nucleon sector, and can be eliminated from consideration by 
what amounts to a different choice of the zero of energy. This is done explicitly by use 
of a field redefinition. As first pointed out by Jenkins and Manohar p!0], the field 
can be re-expressed in terms of a "heavy-baryon" field ipHB^ that has the straight-line 
propagation of a nucleon of mass M removed from its dynamics: 

ip^x) = iPnr{x) exp{—iMv • x). (8) 

If the nucleon is at rest then the choice v = (1,0) is appropriate, although physical 
results should be independent of the choice of v. We also define the analogue of the 
non-relativistic spin 

i 

= (9) 

If we choose v = (1,0) and consider the case of four space-time dimensions then the 
spatial components of S obey the algebra: 

{s:s^} = ls,; [5^S^1=«,.S^ dO) 

Thus a representation for them is provided by the assignment S'^ = SV^? with the 
block-diagonal Dirac matrix constructed out of the Pauli matrices a\ 

Here we focus on the leading pieces of the heavy-baryon Lagrangian, i.e. the 
dominant terms that emerge from Eq. ([g]) in a p/M expansion (with p the three- 
momentum of the nucleon state). The key observation is that the operator structures 
in Eq. ^ do not couple the upper and lower components of t/;^ to one another. Hence, 
the leading behavior of all Dirac bilinears in the p/M expansion can be obtained by 
assuming that ips = H with H obeying ylH = i/, i.e. assuming ips contains only two 
degrees of freedom. We identify these as the degrees of freedom associated with the 
usual non-relativistic spin operator. 

When written in terms of bilinears of i7, the operator structures in Eq. ^ involve 
only the four- vectors v and since: 

Hj^H =[v^]HH (11) 
H^^^.H = 2HS^H. (12) 

These identities, when applied to the version of Eq. ^ with Eq. ([s]) inserted into it, 
produce the standard leading-order heavy-baryon Lagrangian: 

= H [t{v ■ D) + QAU^S^] H. (13) 

The structure of the propagator for the H field can be understood by writing the 
%l) field's total four-momentum, diS k = Mv +p. Then, in terms of the small residual 
("off-shell") momentum the Dirac propagator may be written: 

^ - M ^ 2Mv -p + p^ V -p 2 "" ^ ' 
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The first term is this expansion is then the propagator for the H field. In a typical xPT 
loop graph for a single-baryon process the baryon emits and absorbs a pion, and so it 
is easy to see that the loop graph will be dominated by contributions from momenta 



for which all four components of p are of order 171^^;. The . . . terms in Eq. (14) are 
then suppressed by m^/M. Exceptions to this rule occur in the vicinity of kinematic 
thresholds, but, as long as such thresholds are not nearby, the leading effects of the 
loop graphs in xPT can be obtained by replacing the full baryon propagator by the 



heavy-baryon one (14). 



The corrections of higher order in 1/M can be systematically derived by writing 
ijB = H + h (15) 

with 

= H] ]^h = -h, (16) 

and then calculating the Lagrangian for the field H that results from eliminating the 
degrees of freedom associated with h. This treatment of all 1/M corrections to the 



leading-order heavy-baryon Lagrangian (13) has been carried out by Bernard, Kaiser, 



Kambor, and MeiBner to 0{P^) [llj, Fettes, MeiBner and Steininger to O(P^) f27j, and 
Fettes, MeiBner, Mojzis, and Steininger to O(P^) [28j. At fourth order a large number 
of operators in the non-relativistic theory of which ([l3]) is the leading Lagrangian have 
coefficients that are fixed by the requirement that S-matrix elements at low energies 
agree with the S-matrix elements of the relativistic theory. 

Such considerations regarding the low-energy consequences of Lorentz invariance 
generate the first three terms in the second-order Lagrangian (for further details the 
reader is referred to Ref. [14j). But, since this is an effective theory, the second-order 
Lagrangian contains all operators that are consistent with the symmetries and have 
two derivatives or one power of the quark-mass matrix (since rriq ~ and p ~ ttIt^). 
There are six such independent operators in the isospin-symmetric limit rriu = rud- The 
complete second-order Lagrangian is then [14j: 

= H{^[{v.Dr-D.D)]- '^{S .D,v.u} + c,Tvx+ 
V -uf + C3U-U + (c4 + [5'^, S^]u^u'' 




} H, (17) 



with 



C = eF^,{u^Zu + uZu^), (18) 

the electromagnetic-field-strength tensor (containing both isovector and isoscalar parts) 
that transforms in the appropriate way under the chiral group. In this formulation of 
hadron dynamics the values of two of the low-energy constants (LECs) that appear in 



(17) have already been set to the measured isoscalar and isovector nucleon magnetic 
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moments, tvg and tvy [H The other LECs, C1-C4, encode the effects of more massive 
baryon and meson resonances (A (1232), Roper, p-meson, . . . ) on low-energy dynamics. 
The terms hsted in Eqs. Q, (13), and (17) are aU that is required for computation 



of complete one-loop effects in HB^PT. These loops are divergent, but the degree of 
divergence is mandated by the power of momentum carried by the vertices, and the 
behavior of the heavy-baryon and meson propagators. Dimensional analysis, combined 
with the usual topological identities for Feynman graphs, shows that a xPT graph with 
L loops, N^^ vertices from the nn Lagrangian of dimension rf, and N^^ vertices from 
the ttN Lagrangian of dimension d has a total degree of divergence: 

D = 2L + 2-A + ^{d- 2)Nr + " 1)^^' (19) 

d d 

with A the number of nucleons. 

This calculation of the degree of divergence of the graph is important, because, 
after renormalization is performed, the finite part of the graph must scale as (p, ttIt,-)^. 
Therefore, provided p ^ ttIt^^ the computation of D tells us the order in xPT at which 
the graph contributes. If renormalization is performed using dimensional regularization 
and a mass-independent subtraction scheme then the effects of any particular HB^PT 
graph occur at one and only one chiral order D. This is an advantage of the heavy- 
baryon formulation: each graph scales precisely as P^, and so has a definite order in 



the chiral counting. Specifically, Eq. (19) says that one-loop graphs with vertices from 
C^^l^ give pieces of the quantum-mechanical amplitude with = 3, while graphs with 

(2) 

one insertion from give effects with = 4. Furthermore, since C^TXTx starts at di — 2, 
and Ct^^ at = 1, the lowest order at which a two-loop graph can contribute is = 5. 
Thus one-loop graphs with vertices from l^^^ and either zero or one insertion (s) from 
— together with appropriate tree-level contributions — give the complete xPT result 
up to 0(P4). 

Of course, the divergent loops at O(P^) and 0[P^) need to be renormalized, and 
this is done with counterterms that occur in d>^]^ and Ci"%. Much effort has gone into 
finding a complete, minimal Lagrangian at these orders (see, e.g. Refs. (23, |^J). I 
will not attempt to do justice to that effort here, since only a few of the many contact 
interactions that occur in >C|f^ and C^^\ enter the results presented below. I want to 
stress, though, that it is through such counterterms that hadronic states with excitation 
energies > m^^ affect xPT's predictions. These states therefore produce effects that are 
slowly varying with respect to the external kinematic parameter |q| or while the rapid 
variation of measurable quantities with |q| and uj is provided by pion-cloud physics. This 
interplay of the chiral physics with the higher-energy hadron dynamics is captured in 
xPT, without the need for any assumptions about the details of that dynamics. 

:|: Strictly speaking, it is the chiral-limit magnetic moments that enter here, but the difference is a 
fourth-order effect. 
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3. Chiral perturbation theory in multi-nucleon systems 

The previous section discussed chiral perturbation theory for the single-nucleon system. 
In subsequent sections we wiU examine the predictions of this theory for photon and 
electron scattering on the neutron. However, "experimental" data on these reactions can 
only be obtained using nuclear targets. In analyzing these data it is beneficial to treat the 
interactions that bind the nucleus within the same framework as that used to describe 
the chiral structure of the nucleon itself. This can be done in chiral perturbation theory 
which not only provides a means to compute the low-energy interaction of photons and 
pions with nucleons, but also is used to calculate the forces between nucleons. This is 
the framework that will be used, e.g. in Section [5| to compute photon interactions with 
multi-nucleon systems, and we provide a brief summary of it here. 

The first attempt in this direction was due to Weinberg [ESI El] • Weinberg pointed 
out that HBxPT could not be employed directly to compute the nucleon-nucleon 
amplitude, because a nucleon in a typical nuclear bound state has 

po - pVM; IpI - m,. (20) 

Thus not all components of the (residual) nucleon four-momentum p are order ttZt,-, in 



contrast to the case inside a typical HB^PT loop graph. Because of the scaling (20) 



a diflFerent organization of the Lagrangian (13) plus (17) is required when treating NN 



interactions. The single-nucleon propagator becomes: 

S{p) = '-r^^ '^"^ (21) 

^ ^ 2M ^ 2M 

i.e. the usual non-relativistic propagator for a free particle of mass M. Weinberg further 



observed that in the regime (20) this non-relativistic propagator scales as M/rri^^ which 



invalidates the degree-of-divergence formula (19). 

However, if we consider only diagrams for A^A^ NN that do not contain an 
intermediate two-nucleon state — the so-called 'two-particle irreducible' (2PI) graphs — 



then all single-nucleon propagators have p^ ^ m^^ and Eq. (19) can be used. The goal 
thus becoomes to compute the 2PI NN amplitude up to a fixed chiral dimension D. In 
order to do this we must supplement C^n and C^^^ by a nucleon-nucleon Lagrangian that 
contains the effects of short-distance (r <C physics in the two-nucleon system. 

This physics appears as a string of NN contact interactions of increasing mass dimension, 
each of which has an LEG associated with it. We then identify the NN potential V at 
order D as the sum of all 2PI A^A^ graphs — loops plus contacts — up to that order. 
(Technically an extra step is desirable here, with the 2PI amplitude transformed into 
an energy-independent potential using e.g. the Okubo formalism [30j.) 
The full NN amplitude is then reconstructed from its 2PI parts via: 

r(p',p;^) = l^(p',p) + / g^y(p-,p") ^,_p...^^ T(p-,p;^), (22) 

where an integral over po has been performed to combine the two single-nucleon 



propagators into the propagator of an NN state of energy E. Since Eq. (22) is just 
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the Lippmann-Schwinger equation, which is equivalent to the Schrodinger equation, 
we now have a non-relativist ic quantum mechanics with the NN potential V. The 
resulting approach has been dubbed "chiral effective theory" (xET) since it is a quantum 
mechanics at fixed particle number, rather than a quantum field theory. Within xET not 
only y, but also quantum-mechanical operators for the interaction of external probes 
(pions, electrons, photons, . . . ) can be computed up to a given chiral dimension D. 
If this is done to the same order for current operators and for V then the relevant 
Ward identities are automatically satisfied. But, more generally, this xPT organization 
of quantum-mechanical operators is a systematically improvable way to compute few- 
nucleon structure and reactions at energies of order 171^^. 

In Ref. [31 J the potential V was computed at the one-loop level using C^^^^ and 
up to one insertion from C^^j^ — reproducing the earlier result of Ref. [32j . The formula 



(19) tells us that when this is done we have the complete result for V up to O(P^). 
Representative diagrams at 0{P^) and 0{P'^) are shown in Fig. [l] In Ref. [31j the 
potential was then used (in Born approximation) to obtain phase shifts in the higher 
partial waves of NN scattering. Because of the centrifugal barrier, these partial waves 
are sensitive mainly to the long-distance part of the force. A good description of a 
number of partial waves with / > 2 was found, with the relative contributions of the 
different chiral orders roughly in agreement with predictions of the chiral expansion. 
This suggests that xPT provides a useful organizing principle for the r ~ 1 / part of 
the NN potential. 

\ ^ ^ 

I N / 

' \ / 

I V 



Figure 1. Two diagrams that contribute to the long-range part of the NN potential 
V. The diagram on the left is 0{P^)^ and generates the usual one-pion exchange 
potential, which occurs in the leading-order piece of V. The one on the right is 0{P'^)^ 
and forms part of the NLO correction to the potential. 

This is a key finding, since it opens the way for a treatment of nuclear physics with 
quantifiable uncertainties. If phase shifts can be obtained up to a fixed order n in the 
XPT expansion parameter P, then that suggests that the remaining, omitted physics 
can impact the phase shift only at a level ^ p^+i g^^t the crucial question of how 
to deal with low partial waves, where the NN interaction is non-perturbative, and in 
principle the nucleons probe arbitrarily short distances, was not tackled in Ref. [31j. 

In the seminal paper of Ordonez et al [32j the NN potential was computed up 
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to 0{P^) and Eq. (22) solved for phase shifts in a number of NN partial waves. (See 
also the improved 0{P^) calculation of Ref. [33].) Such a non-perturbative treatment 
of potentials derived from xPT has also been advocated in Ref. [34l [35l [36l |37j. It 
has the advantage that the long-distance parts of and hence of T, are consistent 
with the pattern of chiral-symmetry breaking in QCD. It also means that for r ^ I/tUt^ 
(= |p| ^ ttZtt-) the xPT counting provides a hierarchy of mechanisms within V: the 
leading-order one-pion exchange is more important than the 0{P'^) two-pion exchange 
which in turn is more important than three-pion exchange, and so forth. 

In practice, in order to solve Eq. (22) with the potentials derived from xPT one 
must introduce a cutoff. A, on the intermediate states in the integral equation, because 
the xPT potentials grow with the momenta p and p^ The contact interactions in 
Cnn should then absorb the dependence of the effective theory's predictions on A in 
all low-energy observables, or otherwise we conclude that xET is unable to give reliable 
predictions. In Refs. [32l [33l [38l [39l HQ] V was computed to a fixed order, and then the 
NN LECs that appear in V were fitted to NN data for a range of cutoffs between 500 
and 800 MeV. The resulting predictions — especially the ones obtained with the O(P^) 
potential derived in Refs. [39l[40j — contain very little residual cutoff dependence in this 
range of A's, and describe NN data with considerable accuracy. 

However, several recent papers showed that xET — at least as presently 
formulated — does not yield stable predictions once cutoffs larger than nip 800 MeV) 
are considered [4ll|42l[43l|44l[45l[^ They argue that the theory described above 

is not properly renormalized, i.e. the impact of short-distance physics on the results is 
not under control. Thus the xET described here cannot be used over a wide range of 
cutoffs. Meanwhile, in Ref. [49j it was argued that the cutoff A should be varied only 
in the vicinity of rrip in order to obtain an indication of the impact of omitted short- 
distance physics on the theory's predictions. Since the short-distance physics of the 
effective theory for p ^ rrip is (presumably) completely different to the short-distance 
physics of QCD itself, considering larger cutoffs that allow momentum integrals to probe 
p > rrip does not yield any additional information regarding the true impact of short- 
distance physics on observables. Furthermore, Ref. [49j followed Ref. [34j in arguing 
that working with momenta in the LSE which are demonstrably within the domain of 
validity of xPT has the advantage that the xPT counting must then apply for all parts 
of V: including the contact operators. This justifies xET as a systematic theory, but at 
the cost of limiting its use to cutoffs <C A <C A^sb- Calculations of NN scattering 
using A's in this window have been performed in Refs. j32[ [33l |38l |39l HQ] with, as noted 
above, considerable phenomenological success. 

Consistent three-nucleon and four-nucleon forces have also been derived in xET [501 
ISTl [52] . The three-nucleon force enters first at O(P^), and the four-nucleon force at 
O(P^). This provides an EFT-based explanation for the observed hierarchy of nuclear 
forces: two-nucleon forces provide the bulk of nuclear properties, with three-nucleon 
(3N) forces accounting for small but crucial corrections, and four-nucleon (4N) forces 
almost negligible. 3N and 4N forces have been implemented in calculations of the 
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"low-cutofF' type described in the previous paragraph. For instance, the 3N force 
of Ref. [51j, in conjunction with the two-nucleon force of Ref. [39j, yields a good 
description of nuchdes up to ^ = 13 [53j. An 0{P^) calculation also explains a variety of 
neutron-deuteron scattering data for neutron beam energies up to about 100 MeV [51j. 
Similar success have been achieved for electron-deuteron scattering [54l [55l [56] , Helium- 
4 photoabsorption [57], neutrino-deuteron scattering and breakup [58], pion-deuteron 
scattering at threshold j25l [53 j^, pion photoproduction [62], and — as will be 
discussed extensively in Section [5] — coherent Compton scattering from the NN system. 
From all of this evidence we conclude that structure and reactions in few-nucleon systems 
can be well understood using xET — provided the cutoff is kept below about 800 MeV. 



4. Neutron electromagnetic form factors 

The neutron, as its name implies, is an object with zero nett charge. This does 
not, though, preclude the existence of charged substructure within the neutron. The 
distribution — and to some extent, motion — of charges within the neutron is encoded in 
its form factors. Since the neutron is a spin-half particle there are two of these: electric 
and magnetic in a non-relativistic formulation. These form factors carry the imprint of 
the chiral dynamics of QCD through its impact on the neutron's internal structure. 

The chiral structure of the neutron would be revealed in electron scattering from 
neutron targets — if such experiments were feasible. The gedankenscattering reaction is 
pictured in Fig. [2} where the one-photon-exchange approximation has been assumed. 
The amplitude is then written as 

= (23) 

where q = k — is the four-momentum of the exchanged photon (note < for 
electron scattering) and jj^^ is the leptonic current: 

j;r^ = ti(k')7M«(k). (24) 
The most general form of the neutron relativistic current is the usual: 



«(p). (25) 



2M 

While in the non-relativistic framework being employed here we have: 

j; = d'. {gW)v, + l^GUq')[S,, S,]q^] U, (26) 



where is, in practice, a Pauli spinor [§| and and are the neutron electric and 
magnetic form factors. 

The following discussion is based on the HBxPT analysis of Refs. [HI [63]. For a 
parallel discussion within the context of dispersion relations and models on the role of 
the nucleon's pion cloud in determining electromagnetic form factors see Refs. 



It is a four-component spinor whose lower comonents are zero for the standard choice v = (1,0) 
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k' k 




P' P 



Figure 2. Scattering kinematics for electron- nucleon scattering, in the one-photon 
exchange approximation. In this diagram the nucleon is denoted by the triple line, 
and the electron by the single line at the top of the figure. The blob represents the 
photon- nucleon interaction, ieJ^. 



^EM \ ^V. 



^EM 



For our purposes it is most straightforward to begin with the isoscalar and isovector 
form factors. These are related to the proton and neutron form factors by: 

^E,M^ (27) 

'E,M- (28) 

This decomposition is useful because the leading pion-loop effects occur only in the 
isovector form factors. The quantum numbers of the pion guarantee that the first 
contribution of the pion cloud to the isoscalar electromagnetic form factors occurs only 
at the two-loop level. The one-loop calculation was first carried out in Ref. [llj, and 
was extended and refined in Ref. [63|. Combining the 0{P^) diagrams shown in Fig. 
with the relevant tree graphs yields the expression 



Gl = l 



(47rA)2 



M 



-A 



2 1 

L3^^ + 3J 



(47r/, 



U? Jo 



dx 



{3g\ + l)ml - q'x{l - x){5g\ + 1) log 



log 
■ q'^x{l 



X) 



mt 



(29) 



Here we have worked in the Breit frame, and so = — |qp = Q^, with q = — p the 
three-momentum transfer to the struck neutron. 

The first term, "1" is guaranteed by charge conservation. The second piece gives 
the isovector charge radius, according to the usual definition: 



1 



\{{rlf)q' + 0{q'). 



(30) 



The short-distance contribution to the isovector charge radius, , renormalizes the 
divergent loop integrals. It arises from a term in C\l^ of the form 



B.oHD-f^XH. 



(31) 
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Figure 3. Loop diagrams to be computed to obtain the O(P^) result for the single- 
nucleon form factors in HBxPT. Ah vertices are from C^^j^ and C^}. Picture adapted 
from Ref. |63|. Used with permission of the authors and Elsevier. 



The formula (29) is obtained using dimensional regularization with the modified 
minimal subtraction of xPT for the calculation of the loop integrals. The physical 
result for G^^ is not dependent on the choice of regularization and renormalization 
scheme, but the finite parts of Bi^^jj) will in general be so. The ji dependence is not, 
however, affected by the particular variant of minimal subtraction employed, since the 
II dependence of Bi^^ji) must cancel the ji dependence of the logarithimc term in the 
curly brackets that contributes to the isovector charge radius. 

Note that Bi^ only affects the charge radius, and all higher moments of the charge 
distribution are given — at 0{P^) accuracy — by the physics of the pion-loop graphs 
shown in Fig. [s] The presence of this short-distance contribution to ((r^)^) means 
that the isovector charge radius cannot be predicted in xPT. In fact, the corresponding 
calculation for the isoscalar form factor gives purely short-distance effects: 

2 

I -4Si 



M 



(32) 



Here Bi is the coefficient of the ^NN contact interaction that is the isoscalar analog of 



Eq. (31). However, since there is no loop effect at this order Bi is a finite LEG. 

Thus the prediction for the neutron electric form factor, up to O(P^), in HB^PT 

is: 



1 f.n 

(47rA)2 ' 



1 



4M2 2(47r/^)2 JO 
q^x{l - x){bg\ + 1)1 log 



dx 



{?,gl + l)ml 

— q^x{l — x) 
ml 



(33) 



with B^ a cobination of 5i, 5io, and finite terms with fixed coefficients. The measured 
value of the neutron charge radius squared is ((r^)^) = —0.115(4) fm^ |ESl l67] . 



Comparing this with the result of Eq. (33) reveals that, numerically, the second term 
on the the first line of Eq. (33) — the Foldy term — dominates the overall result, and B^ 
is unnaturally small. The non- analytic behaviour beyond 0{q^) encoded in the integral 



in Eq. (33) is a prediction of xPT at this order, and arises because of the spontaneously 



broken chiral symmetry in QCD. 
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In the case of the magnetic form factor the only counterterm necessary is the one 
already written above in C)^^. It renormalizes the anomalous magnetic moment, i.e. 
the value of the form factor at = 0. Consequently, for the neutron, we find: 

= + 9\ ^^^y \l [jml-q^x{\-x) - . (34) 

Expanding in powers of we find that the magnetic radius of the neutron blows up 
in the chiral limit as m~^/^ . This is to be compared to the log(mg) dependence of the 
electric radius seen in Eq. (29). These chiral-limit divergences of observables — together 
with the associated pre-factors — are model-independent predictions of xPT, since they 
result from the long-distance pieces of the relevant quantum loops. 

In Eq. (34) short-distance physics contributes only to n^. All other moments of 
the magnetization distribution are dominated by pion-cloud physics, and have only 
small {0{P'^) and beyond) contributions from other baryons and mesons. Of course, 
this statement applies only to the regime |q| ^ m^^. Once |q| 600 MeV a variety 
of other dynamics needs to be taken into consideration, since at that momentum 
scale the electromagnetic probe can excite shorter-degrees of freedom in the neutron's 
electromagnetic structure. 

The 0{P^) HBxPT predictions for and are plotted in Fig. [4| together with 
some "data" for both — on which, more below. Also shown is a recent dispersion-relation 
fit to a large database of electron-scattering experiments [67j . 

The upper panel of Fig. [i] shows that the 0{P^) result for has the same feature 
as this dispersion-relation description, namely a rise at low-Q^, followed by a turn 
over. But the HBxPT prediction turns over too quickly. However, adjusting the input 
neutron charge radius to a slightly larger value produces a marked effect. And since 
{{^eY) results from a sizeable cancellation between isovector and isoscalar pieces of 
short-distance physics, this finding suggests that the neutron electric form-factor data 
can be accommodated within a "natural" scenario for that higher-energy dynamics. 
Indeed, the difficulty of describing G^ accurately within HB^PT stems partly from these 
cancellations, which lead to a G^ that is small throughout the entire realm of xPT's 
applicability. This means that higher-order contributions will have a disproportionately 
large impact on the final answer. 

The neutron magnetic form factor, which is order one in this domain, is more fertile 
ground for xPT. And indeed, the 0{P^) prediction for G^ has the right general shape, 
but yields a magnetic radius of 

which is a little smaller than the magnetic radius obtained from the dispersion-relation 
analysis. However, it must be remembered that there is an 0{P^) short-distance 
contribution to this quantity. Including that contribution in the calculation of G^ yields 
the dotted curve shown in Fig. [ij This can be considered a partial O(P^) calculation 
of G^^ since it does not include the additional loop mechanisms that are present at 




Figure 4. Results for (upper panel) and (lower panel) obtained in HBxPT 
up to O(P^) (green dot-dashed line) as compared to the dispersion-relation fit of 
Belushkin, Hammer, and MeiBner (red line, with error bands represented by the shaded 
region) [67]. The (green) dotted line in the upper panel corresponds to an O(P^) 
calculation, but with an input value of {(r^Y) = —0.155 fm^. The measurements 
represented by the circles are from Refs. [68l |69], while the squares are the recent 
BLAST results [70 . The diamonds are from the analysis of data on the deuteron 
quadrupole form factor of Ref. |71J . In the lower panel the data included in the fit of 
Ref. [67^ is shown. The circles are results from deuteron breakup measurements (for 
details see Ref. [64]). The triangles are from the ^He(e, e') measurement of Ref. [72] . 
The dotted line in the lower panel is a partial O(P^) HBxPT calculation of in 
which a 7*nn contact interaction is included so as to bring the result for the neutron's 
magnetic radius into agreement with the central value obtained in Ref. [67j . 

this order. But even this partial 0{P'^) computation shows that the shift due to short- 
distance effects is about 20% at = 0.1 GeV^ — a magnitude that is consonant with this 
being a contribution of 0((3^/A^gg) to G^/tvn- In fact, the O(P^) contact interaction 
that shifts the magnetic radius to the empirical value is associated with a momentum 
scale of order 1 GeV. The LEG that is the coefficient of this contact interaction is 
certainly natural with respect to A^sb- But by ^ 0.25 GeV^ the O(P^) contribution 
is as large as the 0{P^) result — a clear sign that the xPT expansion is breaking down. 
The results for therefore confirm the xPT picture that the dominant effect in the 
neutron's magnetization distribution is the pion cloud, while also showing that this 
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picture is limited to momentum transfers below 500 MeV. 

The results reported here are rigorous within HBxPT. Ref. [63] showed that the 
inclusion of Delta(1232) degrees of freedom in the effective field theory does not improve 
the agreement with data — the short-distance physics in these observables is associated 
with higher-mass states. The impact of corrections associated with nucleon recoil on 
neutron form-factor predictions is discussed in Refs. [731 El!- These corrections are 
appreciable in the case of G^, and there the q/M series can be reorganized into a 
form that improves the agreement of the xPT predictions with data. The tree-level 
contribution of vector mesons is also considered in Refs. [731 El] and it assists in the 
description of G^. 

The "data" points represented by circles and squares in Fig. |4] are obtained via 
electron-induced breakup reactions on deuterium targets. It needs to be emphasized 
that xPT calculation of these reactions would allow a completely consistent treatment 
of the pion dynamics that generates the higher-order terms in the neutron form factor 
and the pion-exchange currents which play, e.g. an important role in 7^ np at 
threshold ^5]. These exchange currents occur at 0{P^) in the irreducible kernel for 
the interaction of real and virtual photons with the NN system, and so have the same 
chiral order as the single- nucleon loop effects that have been highlighted in this Section. 
But, no comprehensive analysis of ed e'np has been performed within xET. Such an 
analysis, and its extension to the three-body system, would allow a consistent treatment 
of the chiral structure of the neutron and the target nucleus. It would also facilitate 
reliable estimates of the theoretical uncertainties in the neutron form-factor extraction. 

The diamonds in the upper panel of Fig. [4] are obtained from an analysis of the 
deuteron's quadrupole form factor, Gq In this case the AV18 potential [76j was 

used to generate the deuteron's wave function, and the important correction to the 
deuteron's charge operator [55l [77l [78] was included in the calculation. However, this 
calculation fails to reproduce the deuteron's quadrupole moment, and this could have 
some consequences for the Q^-dependence of the deuteron Gq [56j. xET describes 
data on the ratio of deuteron charge to quadrupole form factors Gc/Gq well up to 
momentum transfers of about 0.3 GeV^ [Mj- A xET extraction of from Gg, and 
possibly from the A{Q'^) data of Ref. [79] and subsequent studies, could perhaps provide 
useful information on G^ in this \ow-Q^ regime. 

5. Compton scattering from the neutron 

The neutron has zero charge, but its composite structure means that it has a non-zero 
Compton scattering amplitude. The standard decomposition for the neutron Compton 
amplitude involves six non-relativist ic invariants: 

T^N = Me • e + A2e-ke-k' + iA^a • {e x e) + iA^a • {t x k) e - e 
+ lA^G • [{e xk)e't -{ex k') e 'k]+ iAqg • [{e x t) e - t - {e x k) e -k], 

(36) 
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where k and e (k^ and e^) are the three-momentum and polarization vector of the 
incoming (outgoing) photon. The A^'s, i = 1 . . . 6, are scalar functions of photon energy 
and scattering angle. 

At photon energies ^ m^j^ Compton scattering from the neutron is driven by photon 
interactions with the neutron's pion cloud. The dominant effects come from the 0{P^) 
mechanisms depicted in the diagrams (a)-(d) of Fig. [5] and the first line of Fig. |6] These 
loop diagrams individually contain divergences, but the divergences cancel in the sum. 
Indeed, they must do so, because no gauge-invariant jn contact interaction can be 
constructed at 0{P^): at this level of accuracy there is no short-distance contribution 
to neutron Compton scattering. Consequently, up to 0{P^) there is a xPT prediction 
for this process. It, together with the associated proton result, is given by the following 
(Breit-frame) expressions for Ai-Aq 
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with 

S = yjl - T^z^ - t{l - zyx{l - x), R = ^Jl-t{l- zyx{l-x). (39) 

Here and below formulae are given in Coulomb gauge, for real photons, so the 
polarization four- vectors and photon four-momenta satisfy: 

P = 0; = 0; V'e = 0] v-e = 0. (40) 




(d) (e) (f) 



Figure 5. Tree diagrams that contribute to Compton scattering in the e-v = gauge 
up to O(P^). Small dots are vertices from >C^^^, larger dots are vertices from >C^^^, the 
sliced dots are vertices from >C^^^ and the hatched dot is a vertex from >C^^^. Crossed 
graphs and graphs which differ only in the ordering of vertices are included in the 
calculation, but are not shown here. Figure from Ref. [82 , reprinted with permission 
of Elsevier. 

Note the impact of the Goldstone bosons of QCD — the pions — on the neutron 
Compton amplitude. The theory predicts that the Compton amplitude contains a cusp 
ai uj = m^. That cusp is a consequence of the opening of the photoproduction channel 
at that energy [Hj but its consequences are felt even for cu well below rrij^. Thus Compton 
scattering — alone among reactions with real photons — probes chiral dynamics at photon 
energies well below the pion threshold. 

The most famous manifestation of this chiral dynamics in neutron Compton 
scattering occurs in the neutron polarizabilities. The low-energy expansion of the Ai^s 
defines the electric and magnetic polarizabilities a and f3 as well as the four "spin 
polarizabilities" 71^74. The former occur as coefficients of the cj^ terms in Ai and ^42, 
while the latter are coefficients of the terms in A^-Aq — once the Born terms are 
accounted for. For the Breit-frame amplitude such an cj-expansion gives (keeping terms 
up to 0(1/ M^)) [HUES]: 

M^.0)=-^ + ^^ ((Z + + cos e)-Z'){l- cos 9) 

+ 47r(a + /3 cose)6j2 + 0(cj^), 
A2{uj, 9) = -^^{^Z + /^)6j2 COS 9 - iirPuj^ + 0(6J^), 

II In HBxPT the threshold is not at precisely the right photon energy, but this can be corrected through 
a resummation of higher-order terms that is motivated by a modified power counting in the vicinity of 
a threshold [80|. 




Figure 6. Diagrams which contribute to nucleon Compton scattering in the e ■ v = 
gauge at 3rd (i-iv) and 4th (a-s) order. Vertices are labeled as in Fig. [S] Fi gure 
reprinted from Ref. [82], with permission from Elsevier. 



^,{2{Z + 2k) -{Z + Kf cos 6) + Af 

+ 47rcu^(7i - (72 + 274) cos 6') + 0{u^), 
{Z + Kf + Af + 4716^^72 + 0(a;5), 



2M2 



A^{u,d) = 



2 

2M2 



+ + +47ra;373 + 0(6^5), 



(41) 



with ^3°-y4.Q° the contributions from the pion-pole graph, Fig. [5j 

At O(P^) the xPT predictions for these polarizabihties are straightforwardly 
obtained from a Taylor series of Eq. (37) in powers of uj/rriT^. The following results 
for the polarizabihties can be obtained ^SlIlTj by matching Eq. (37) to Eq. (41): 
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- 1.1 X 10"Mm^ 



3847r3/2m2 - - ' (42) 

As with the charge and magnetic radii electromagnetic polarizabihties diverge in the 
chiral hmit of QCD — with the coefficient of the divergence predicted by HB^PT. That 
divergence is particularly dramatic l/rriq) in the case of the spin polarizabihties, and 
this could have interesting consequences for lattice simulations of these quantities [84j . 
And the HB^PT result for the 7-nucleon amplitude actually gives much information 



besides the predictions (42), since HB^PT predicts the full dependence on the parameter 

But and can be extracted from low-photon-energy proton Compton-scattering 
experiments. Such measurements are particularly revealing because at O(P^) a and 
f3 receive their first correction from short-distance physics. They thus result from 
a fascinating interplay of the pion-cloud dynamics that yields the dominant effects 



of Eq. (42), and shorter-distance, higher-^PT-order mechanisms that correct those 



predictions. 

Elsewise at O(P^), nucleon-pole graphs [5](e) and the fixed-coefficient piece of the 
seagull depicted in Fig. [5](f) give further terms in the 1/M-expansion of the relativistic 
Born amplitude. There are also a number of fourth-order one-pion loop graphs (graphs 
a-r of Fig. p). All of the graphs pa~l6r involve vertices from C^^^^. In particular, the 



from C^^l 



expressions for diagrams [6|i-[6|: contain the LECs Ci, C2, and C3, and in diagrams |6f 
the nucleon anomalous magnetic moments enter. 

Divergences in graphs [6^-[6]f, [6|i-[6|^, [6|), and are renormalized by counterterms 
The overall result — first derived in Ref. jSOj — is 




12/X5)r3 



(43) 



+ 51 + 24/x,r3 



+ 2c3 - 4ci 



8M^ 

+ 5/?^(/x)+r35/?^(/x). (44) 

The dependence on renormalization scale, /x, in the isoscalar and isovector LECs 
5a^'^ and 5f3^'^ cancels that of the loops. These LECs encode contributions to the 
polarizabihties from mechanisms other than soft-pion loops, i.e. short-distance effects. 
In contrast, the spin polarizabihties, as well as all higher terms in the Taylor-series 
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expansion of the Compton amplitude, are still predictive in xPT at 0{P^). (The full 
amplitudes up to 0{P'^) are not written here; they can be constructed from the details 
given in the Appendices of Ref. [82j.) 

Use of the O(P^) amplitude allows a good description of 7p data in the kinematic 
range cj, \f\t\ < 180 MeV (see Fig. 7). This calculation includes the effects (37), and 

The 



also the additional 0{P^) effects in Fig. [5] and |6] that were computed in Ref. [81 
only free parameters are and (3^. A fit to data in this kinematic range yields [82^ 185]: 

= (12.1 ±1.1)1^1 X 10-^fm^ 

= (3.4 ± l.l)tol X 10-^ fm^ (45) 

Statistical (1-cr) errors for a two-parameter fit are inside the brackets, while the second 
set of errors are an estimate of the uncertainty due to the effect of higher-order terms. 
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Figure 7. Results of the 0{Q^) HBxPT best fit to the differential cross sections for 
Compton scattering on the proton at various angles, compared to the experimental 
data [86l [87]. The gray region is excluded from the fit. The magenta diamonds are 
Mainz data [87]; the other symbols are explained in Ref. [81]. Figure reprinted from 
Ref. [82 , with permission from Elsevier. 



This success in 7p scattering prompts an important question about neutron 
polarizabilities: how different are they from those of the proton? In particular, the 
leading pion-cloud mechanisms in the first line of Fig. |6] all yield exactly the same 
contributions for the neutron and the proton. So too do the dominant mechanisms 
associated with excitation of the Delta(1232) ^[j^. So, a measurement of a non-zero 
value of the isovector quantities and f3^ (defined as in Eq. (44)) would be indicative 
of other (i.e. non-leading-Delta-excitation) short-distance effects in the polarizabilities. 

In order to obtain information on neutron polarizabilities we must analyze Compton 
scattering from light nuclei. At photon energies ^ 171^^ the approach to reactions 
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described in Sec. [3] tells us that the Compton amplitude for these nuclei is formed 
by computing 

T^A = {^a\T^n + T^nnI^a) (46) 

where I^Pa) is the wave function of the nucleus, and T^n is the single- nucleon Compton 
amplitude. T^nn is the amplitude for 7A^A^ ^NN^ with all NN interactions after 
(before) the departure (arrival) of the outgoing (incoming) photon amputated from it. 
If these amplitudes are computed up to order n in xPT, and the wave function \ipA) 
is obtained from a xPT potential that is also computed at that order, then all the 
consequences of chiral-symmetry breaking for the 7A amplitude have been accounted 
for — up to the 0{P^^^) corrections occurring at higher orders in the xET expansion. 

The xPT amplitude T^nn was first computed in Ref. [89j, and includes the diagrams 
shown in Fig.jsj At 0{P^) it contains the two-nucleon analog of the mechanisms depicted 
in graphs (i)~(iv) of Fig. |6| These isoscalar exchange currents are sizeable at the ^ 100 
MeV photon energies of interest to us here. 




Figure 8. Two-body diagrams for the irreducible ^NN ^yNN kernel, T^atat, at 
O(P^), sandwiched between deuteron wave functions which are represented by the oval 
blobs. Permutations are not shown. Figure reprinted from Ref. [82^, with permission 
from Elsevier. 



The xET computation of 7d scattering has received extensive attention, with an 
0{P^) computation in Ref. |82llH5]. That calculation showed that short-distance eflFects 
in the 7d reaction are not described with enough precision at 0{P^) in xET to allow 
an accurate extraction of polarizabilities from existing experimental data. The short- 
distance physics that is missing is of two types. 

First, the Delta(1232) starts to have an impact beyond its effects on the "static" 
polarizabilities, (44), already at 80 MeV — well below the upper energy of the 
data set. This effect is particularly noticeable at backward angles, where magnetic 
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scattering is known to lead to significant Delta excitation in the 7p case [20l EH [22] . 
Since the value of ^J\t\ ^150 MeV for some existing data, it is no surprise to find 
Delta(1232) effects that are nominally 0{P^) in HBxPT have a ^ 10% effect on the 
7d differential cross section. This mars our attempts to describe 7d data in HB^PT at 
energies of order 100 MeV, and translates into ^ 100% uncertainties in (3^. 

This deficiency was rectified in Ref. [90j, where the leading effects of Delta(1232) 
excitation were added to the 7d calculation. This requires the use of a new effective 
theory, in which the Delta(1232) is included as an explicit low-energy degree of 
freedom [iHl [HI [21] . That new EFT provides an accurate description of higher-energy 
data in both 7N and 7d scattering — as shown in Fig. [9] for the 7d case. 

Second, there are short-distance, two-body Compton mechanisms that are 
suppressed by two powers of P, compared to the dominant 0{P^) piece of T^nn- But, 
when xPT is applied to photon-nucleus scattering, the expansion parameter is the larger 
of }J(jj/M and m^^/ \/McJ. Therefore these higher-order pieces can easily give a 10% effect 
in the amplitude at 100 MeV. This was the size of differences seen in Ref. [82] when 
deuteron wave functions with different implementations of short-distance NN dynamics 
were used in the calculation. 
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Figure 9. Results for Compton scattering from the deuteron at photon lab. energies 
49 MeV and 94.2 MeV |91j. Data shown are from Refs. [92l|93] respectively. The long- 
dashed curve is an O(P^) calculation a la Ref. [89^, while the short-dashed curve is the 
result found in one of the fits of Ref. [82]. The solid line supplements the O(P^) result 
with the dominant effects of Delta(1232) excitation, and corresponds to the central 
values in Eq. (48). Figure courtesy H. GrieBhammer. 



This difficulty can be circumvented by recognizing that in the regime of very 
low photon energies a different power counting is required. The power counting for 
uj ^ m\/M was implemented in Ref. 



In that domain Eq. (46) is modified to: 



\rpirr _|_ rpirr 



+ fGJ + JGf\ijA), 



(47) 



where J is the operator for photo-absorption by the two-nucleon system and G is the 
fully interacting two-nucleon Green's function. In Eq. (47) the operator T^)^ includes 
only the one-nucleon-irreducible part of the 7N amplitudes, and a similar definition 
applies to T^^^n- This avoids double counting with the last two terms on the right-hand 
side of Eq. (47). If J, X!^, T^^nn-) the NN potential V are all computed up to 
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the same order in xPT then the calculation of Eq. (47) is manifestly gauge invariant. 
Consequently, the correct Thomson limit is automatically obtained, and a significant 
reduction in the sensitivty of T^a to details of the short-distance physics in the NN 
system results. 

This repairs the missing physics in the calculation of Ref. f82] because the part of 
J^GJ + JGJ^ that was not included in T^n and T^nn of Eq. (46), while higher-order in 
XET, is not gauge invariant by itself. Once a gauge-invariant calculation is performed the 
short-distance NN dynamics can only impact the amplitude at 0{P^) and beyond. 
The calculations of Ref. [94] show that these effects are very small, with the uncertainty 
in the cross section due to the use of deuteron wave functions \ipd) that contain different 
short-distance physics being < 2% in the experimentally relevant energy range. 

That range presently extends from uu ^ 45 MeV to uu ^ 95 MeV, and over angles 
from about 30^ to 150^. An extraction of neutron polarizabilities that is as accurate as 
that of Eq. (45) for the proton case therefore requires the formalism of Ref. [H^, as well 
as some discussion of Delta(1232) effects. But, it is important to note that for uu = 50- 
100 MeV the dominant effects in this computation arise from excitation of the pion 
cloud: the pion cloud of the nucleon that gives the dominant contribution to the electric 
polarizability, and the pion cloud of the deuteron that produces the exchange currents 
shown in Fig. [8| Ref. [94j analyzed the data on elastic scattering and extracted the 
result [91j: 

a' = (11.3 ± 0.7 ± 0.6 ± 1.0) x 10"^ fm^ 

(3' = (2.8 ± 0.7 T 0.6 ± 1.0) x 10"^ fm^ (48) 

for the isoscalar polarizabilities. (The first error is the 1-a statistical error, the second 
comes from the use of the Baldin sum rule to constrain the fit, and the third is an 
estimate of the potential impact of omitted mechanisms on the result.) The quality of 
the fit is also good, see Fig. [9] for examples at photon laboratory energies of 49 MeV 
and 94.2 MeV. The ongoing experiment at MAX-Lab in Lund, Sweden [95], as well 
as future experiments at the High-Intensity Gamma-ray Source (HI7S) at the Triangle 
Universities Nuclear Laboratory, will provide much-needed additional data to test this 
theoretical description, and improve the extraction of and /?^. 

Recently the first computation of elastic Compton scattering on the three-body 
system was completed j96l |97j . This was also the first consistent xET calculation of 
an electromagnetic process in the three-nucleon system. In this case the focus was 
on photon energies in the 60-120 MeV range, and the formula (46) was used. This, 
then, was a first, exploratory, study that sought to examine what information Compton 
experiments using a Helium-3 target could give us about neutron polarizabilities. Once 
again the T^n and T^nn were computed up to O(P^), and a variety of potentials were 
employed in the calculation of \ipHe^)- 

The first conclusion from this study is that Helium-3 will have a significantly larger 
cross section for coherent Compton scattering than deuterium, due to the presence of 
two protons in the nucleus. Furthermore, the Thomson terms of these two protons can 
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interefere with the neutron polarizabihty pieces of the 7N amphtude, so these pieces 
have a larger effect than in scattering — at least in absolute terms. 

In addition, polarized Helium-3 is unique and interesting because it seems to behave 
as an "effective polarized neutron" in Compton scattering. The Compton response is — 
at least up to 0{P^) — dominated by the nuclear configuration where the two protons 
are paired in a relative S-wave, and the spin of the nucleus is carried by the neutron. 
And two-body currents at this order are largely spin independent. These two facts lead 
to 7^He double-polarization observables that look very similar to the asymmetries that 
would be measured were the same double-polarization experiments done on a neutron 
target. 

Figure [To] shows predictions for A^, which is defined as a difference of differential 



cross sections with circularly polarized photons and the target polarized parallel and 
anti-parallel to the beam direction, as a function of the scattering angle of 120 MeV 
photons. Here different neutron spin polarizabilities are implemented in an otherwise 
strict, O(P^), 7He^ scattering calculation. In each panel one of the four neutron spin 
polarizabilities is varied. This shows that a measurement of at this energy should 
allow extraction of the combination 7]^ — (72 + 274) cos^ (see Ref. [97j for details). 




Figure 10. The four panels above correspond to the O(P^) result for A^, with 
an additional variation of, respectively, Aji^ (top-left panel), Aj2n (top-right panel), 
Ajsn (bottom-left panel) and Aj^ri (bottom-right panel). The calculations are done in 
the cm. frame at 120 MeV. The solid (black) curves correspond to the full unperturbed 
O(P^) results. The long-dashed (blue) curves correspond to Aj^ = — 7j^(0(P^)), dot- 
dashed (red) to A7f = -7^ (0(P^))/2, dotted (magenta) to A-/^ = 7^ (0(p3))/2 and 
dashed (green) to Aj^ = jin{0{P^). Figure reprinted from Ref. [97], with permission 
from Elsevier. 



Next, in Fig. 



11 



we plot Ax vs. the cm. angle at 120 MeV. Again, the different 
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Figure 11. The four panels above display the O(P^) results for A^^ at a cm. photon 
energy of 120 MeV. Also shown are the results of varying one neutron spin polarizability 
per panel. Legend as in Fig. 10 Figure reprinted from Ref. [97 , with permission from 
Elsevier. 



panels each correspond to varying one of the four spin polarizabilities in an otherwise 
strict O(P^) xET calculation. The results suggest that A^^ is sensitive to a different 
combination of neutron spin polarizabilities than is A^. 

Measurement of these obseravbles is a high priority for HI7S, along with polarized 
7P and 7d. The idea is that the jp jp experiments will reveal proton spin 
polarizabilities, which can then be compared to the neutron ones measured in the 7He^ 
experiments. A different combination of neutron spin polarizabilities can be accessed 
through double-polarization measurements on the deuterium nucleus [98j . There is also 
the possibility of coherent Helium-3 Compton scattering measurements at MAX-Lab in 
the near future. This allows us to anticipate the measurement of several neutron spin 
polarizabilities, with the consequent possibility of examining differences between proton 
and neutron 7^'s. There is also progress on the extraction of information on these 
quantities from lattice QCD ^9], and so such experiments could provide important 
insights into the interplay of Goldstone-boson and higher-energy QCD dynamics in 
nucleon Compton scattering. 

6. Conclusion, and some words on other electromagnetic-induced reactions 

Electron scattering from the neutron provides the opportunity to map out its 
electromagnetic structure as a function of Q^, the square of the space-like momentum 
transfer in the reaction. Neutron Compon scattering provides a different picture, as 
we probe neutron structure as a function of photon energy uu. In both cases the rapid 
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variation of response functions with the kinematic parameter is driven by the neutron's 
pion cloud. Slower variation comes from excitations of higher energy. The review of 
Ref. places greater emphasis on the role of those higher-energy degrees of freedom, 
and also covers a wider swath of kinematical territory, than does this article. The 
understanding of the nucleon's electromagnetic response I have presented here is based 
on the chiral symmetry of QCD, and the status of pions as Goldstone bosons of that 
symmetry. However, since chiral symmetry is only an approximate symmetry of the 
QCD Lagrangian, the "slow" variation from more massive degrees of freedom can, on 
occasion, be just as fast as that due to pions, see e.g. the problem of the magnetic 
response of the nucleon in Compton scattering ESj • 

It is critical that we test and elucidate the consequences of QCD associated with 
this physics for both neutrons and protons. Chiral perturbation theory systematically 
implements the interplay of pion-cloud and shorter-distance mechanisms in nucleon 
structure. It predicts definite patterns for proton and neutron observables, e.g. that it 
is the isovector combination of form factors that has the dominant effects due to chiral 
symmetry, or that proton and neutron polarizabilities are equal at leading order. The 
calculations and experiments reported on here show that a picture of neutron structure 
based on chiral symmetry provides a good qualitative understanding of the neutron's 
electromagnetic excitation already at the leading one-loop level. That picture can be 
rendered quantitatively accurate if some care is taken to include the dominant short- 
distance mechanisms, e.g. those due to the Delta(1232), through use of the chiral 
effective field theories developed in Refs. pTOl [TSl [T9l [20] . 

An obvious extension of the tests discussed here is to examine the neutron's response 
as a function of u and simultaneously. This has been done for the proton, through 
(e, 6^7) or "virtual Compton scattering" (VCS) experiments. Chiral perturbation 
theory up to 0{P^) provides a good description of the lowest-Q^ experiment of this 
type At this order, xPT not only correctly predicts the electric and magnetic 

polarizabilities of the proton to within experimental uncertainties, it also predicts the 
spatial distribution of the polarizability response within the proton jl02j . The theory 
also predicts that the neutron's generalized polarizabilities are the same as the proton's 
at this order. Experiments to test this prediction would be very timely, and could be 
accomplished through a VCS experiment on the deuteron, together with a calculation 
of the two-body corrections to 0{P^). More ambitiously, generalized neutron spin 
polarizabilities could perhaps be extracted from VCS experiments of the type recently 
undertaken at Mainz in the proton case |103j . but with the proton replaced by a Helium- 
3 target. 

Of course, the most direct evidence for the chiral structure of the neutron comes 
from pion-production reactions. The amplitude for 777, 7T~p is given to very good 
accuracy by the leading chiral effect: the "KroU-Ruderman" term. In this case 
O(P^), where the leading-loop contributions appear, is two orders beyond leading and 
interesting short-distance effects also enter at this order [104] . The 0{P^) amplitude 
for charged-pion photoproduction on the neutron can be compared to experiment by 
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examining data on pion capture by the proton — data which is in good agreement with 
the xPT result. An approved experiment at the MAX-Lab facihty in Lund, Sweden, wiU 
soon measure the energy dependence of the near-threshold cross section for jd Ti~pp 
to good accuracy [105j . This should provide useful information on the charged-pion 
photoproduction reaction on the neutron, although again, this one-body part of the 
amplitude will have to be disentangled from NN mechanisms. 

Finally, one of the most beautiful pieces of evidence in neutron dynamics regarding 
the pattern of chiral-symmetry breaking occurs in the reaction 772 n^n. This is 
predicted to have a large threshold £"0+ amplitude, in spite of the absence of any neutron 
charge, because of the impact of the two-step process 777, 7T~p n^n — even though 
the threshold for production is below the charged-pion threshold. HB^PT predicts 
the neutron £"0+, based on the one-loop dynamics corresponding to this process, as well 
as the measured proton value, £9+^, and the assumption that the short-distance physics 
in the proton and neutron reactions is the same (apart from obvious isospin factors). 
The result — which includes terms up to 0{P^) in HB^PT — is |l()6j : 

Ef^^ = 2.13 X 10-^m-\ (49) 

HBxPT's prediction for 777. n^n can be tested in n^d. But, a consistent 

calculation of the threshold S-wave amplitude for neutral-pion photoproduction on 
deuterium, E^^ requires the inclusion of the 0{P^) and 0{P^) mechanisms for 7A^A^ 
TT^NN [61[|62]. These diagrams include the double-scattering contributions that are the 
two-body counterpart of the diagrams that in the single-nucleon case lead to the large 
neutron Eq^. xET up to 0{P'^) then predicts: 

Ed = (-2.6 ± 0.2) X + 0.38£J°^. (50) 

The uncertainty comes from the experimental error on Eq^ and the variation in the 
result with different NN-system wave functions. (Note there are other higher-order 
uncertainties that are not included in this error bar.) Using a sequence of measurements 
at energies up to 20 MeV above threshold |107j to obtain E^^ and then using the result 
of Eq. (50) to convert this to a result for the threshold n^n amplitude, we find: 

E^^ = (3.03 ± 0.58) X 10-^m-\ (51) 

with an additional, here unquantified, uncertainty arising from the subtraction of the 
inelastic channel that was employed in Ref. [107J. The reasonable level of agreement 
between the 0{P^) prediction (49) and the measurement (51) shows that we understand 
the chiral dynamics of the neutron, as revealed in ^d ^ n^d at threshold, quite well. 

Unfortunately the situation in ed e'n^d is not quite as clear |108j . For the case of 
threshold pions and low momentum transfer, a successful description could reasonably 
have been expected, since the pion cloud should also dominate this process. However, 
there are puzzling discrepanices between theory and data already for = 0.1 GeV^. It 
is hoped that careful treatment of the NN system dynamics, together with a thorough 
examination of the Delta(1232) mechanisms that enter this process, can explain why 
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the pion cloud dominates vr^d, but is in competition with higher-energy QCD 

dynamics in the case of electro-production. 

Studies that explore these issues will ensure that photo- and electro-production 
continue to provide important windows on the chiral structure of the neutron. A unified 
picture of low-energy and low-momentum-transfer neutron structure as revealed in these 
reactions, as well as neutron Compton scattering (real and virtual), and the neutron 
electromagnetic form factors has been established. This picture — which is systematized 
by chiral perturbation theory — is grounded in QCD and agrees with data in all these 
reactions at the expected level given the xPT order to which calculations have been 
accomplished thus far. Ongoing efforts in both theory and experiment — and continued 
dialogue between them — will elucidate further details of neutron structure, and so 
improve our understanding of QCD in its strongly-interacting regime. 
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